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Abstract. The Wyman's solution depends on two parameters, the mass M and the 
scalar charge a. If one fixes M to a positive value, say Mq, and let cr^ take values along 
the real line it describes three different types of spacetimes. For ct^ > the spacetimes 
are naked singularities, for = one has the Schwarzschild black hole of mass Mq and 
finally for —Mq < < one has wormhole spacetimes. In the present work, we shall 
study qualitative and quantitative features of orbits of massive particles and photons 
moving in the naked singularity and wormhole spacetimes of the Wyman solution. 
These orbits are the timelike geodesies for massive particles and null geodesies for 
photons. Combining the four geodesic equations with an additional equation derived 
from the line element, we obtain an effective potential for the massive particles and a 
different effective potential for the photons. We investigate all possible types of orbits, 
for massive particles and photons, by studying the appropriate effective potential. We 
notice that for certain values of > 0, there is an infinity potential wall that prevents 
both massive particles and photons ever to reach the naked singularity. We notice, 
also, that for certain values of —Mq < cr^ < 0, the potential is finite everywhere, which 
allows massive particles and photons moving from one wormhole asymptotically flat 
region to the other. We also compute the radial timelike and null geodesies for massive 
particles and photons, respectively, moving in the naked singularities and wormholes 
spacetimes. 
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1. Introduction 

The weak equivalence principle of general relativity tell us that massive particles move 
along timelike geodesies and photons move along null geodesies [1]. Despite their 
fundamental importance as one of the principles of general relativity, the geodesies 
also help us learning more about different properties of a given spacetime. A textbook 
example comes from the study of geodesies in Schwarzschild geometry [1]. Without 
actually computing the geodesies, just observing the effective potential diagram, one 
can see that both massive particles and photons can never leave the event horizon once 
they enter that surface. This is the case because the effective potential for both massive 
particles and photons diverges to negative infinity as one approaches the singularity 
located at the coordinate system origin. Therefore, once the massive particles and 
photons enter the event horizon they are accelerated toward the singularity without any 
chance to turn back. Many authors, over the years, have computed the effective potential 
diagram and the geodesies of different spacetimes in order to learn more about their 
properties. In particular, we may mention some important works dealing with different 
black hole spacetimes Two other important gravitational configurations besides 
black holes that may form due to the gravitational collapse are naked singularities 
and wormholes. Some authors have already investigated some of their properties by 
computing effective potential diagrams and geodesies for those space-times [31 HI [5]. 
A well-known spacetime geometry which may describe naked singularities as well as 
wormholes is the Wyman one [6] . Since, to the best of our knowledge, nobody has ever 
studied the Wyman spacetime by means of the computation of its effective potential 
diagrams and geodesies, we decided to do that in the present work. 

The Wyman geometry describes the four dimensional space-time generated by a 
minimally coupled, spherically symmetric, static, massless scalar field and has been 
studied by many authors [6l [3, El El [ini [H]. From a particular case of the general 
Wyman's solution, M. D. Roberts showed how to construct a time dependent solution 
[Hj. The Roberts' solution has an important physical interest because it may represent 
the gravitational collapse of a scalar field. Later, P. Brady and independently Y. Oshiro 
et al. [12] , [13] showed that the Roberts' solution could be derived from the appropriated, 
time-dependent, Einstein-scalar equations by using a continuous self-similarity. They 
also showed that the Roberts' solution exhibits a critical behavior qualitatively identical 
to the one found numerically by M. W. Choptuik [H], studying the same system of 
equations. In fact, the above results confirmed early studies of D. Christodoulou who 
pioneered analytical studies of that model [15j. 

The Wyman's solution is not usually thought to be of great importance for the 
issue of gravitational collapse because it is static and the naked singularities derived 
from it are unstable against spherically symmetric linear perturbations of the system 
[H [To] . On the other hand, as we saw above, from a particular case of the Wyman's 
solution one may derive the Roberts' one which is of great importance for the issue of 
gravitational collapse. Also, it was shown that there are nakedly singular solutions to 
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the static, massive scalar field equations which are stable against spherically symmetric 
linear perturbations (TO]. Therefore, we think it is of great importance to gather as 
much information as we can about the Wyman's solution for they may be helpful for a 
better understand of the scalar field collapse. 

The Wyman's solution depends on two parameters, the mass M and the scalar 
charge a. If one fixes M to a positive value, say Mq, and let take values along the 
real line it describes three different types of spacetimes. For > the spacetimes 
are naked singularities, for cr^ = one has the Schwarzschild black hole of mass Mq 
and finally for — Mq < < one has wormhole spacetimes. Therefore, we have an 
interesting situation where we can study different properties of naked singularities and 
wormholes together through the effective potential diagrams and geodesies of massive 
particles and photons in the Wyman's solution. 

In the next Section, we introduce the Wyman solution and identify the values 
of the parameters M and t] that describes naked singularities, wormholes and the 
Schwarzschild black hole. In Section [3l we combine the geodesic equations to obtain 
the effective potential equation for the case of massive particles. We study the effective 
potential and qualitatively describe the types of orbits for massive particles moving in 
the naked singularities and wormholes spacetimes. In this section, we also compute the 
radial timelike geodesies for particles moving in the naked singularities and wormholes 
spacetimes. In Section HI we combine the geodesic equations to obtain the effective 
potential equation for the case of photons. We study the effective potential and 
qualitatively describe the types of orbits for photons moving in the naked singularities 
and wormholes spacetimes. In this section, we also compute the radial timelike geodesies 
for photons moving in the naked singularities and wormholes spacetimes. Finally, in 
Section [5] we summarize the main points and results of our paper. 

2. Naked singularities, wormholes and the Schwarzschild black hole. 

The Wyman line element and the scalar field expression are given in terms of the 
coordinates {t, r, 9, cf)) by the equation (9) of [8], 



where r varies in the range 2// < r < oo and dfl"^ is the line element of the two- 
dimensional sphere with unitary radius. The scalar field is. 



where a is the scalar charge given by a'^ = rf — and cr^ > — M^. We are working in 
the unit system where G = c = 1. 

In the line element ([1]), the function -R(r), 






(2) 




(3) 
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is the physical radius which gives the circumference and area of the two-spheres present 
in the Wyman geometry. Instead o using the parameters M and a to identify the 
different types of spacetimes described by ([T]), we may also use the parameters M 
and Tj, since they are all related by the equation just below ([2]). Each different type 
of spacetime will have a different behaviour of the function R{r). R{r) has a single 
extreme point which is a minimum located at r = M + r^. Since, the line element ([1]) 
is valid for r > 21], only in the cases where M > t] there will be a minimum inside the 
domain of r. Based on that result we have three different cases for a positive M. As 
a matter of simplicity we shall use, also, the positive parameter A = M/rj in order to 
identify the three different cases. 
Case M<r] or 0<X<1. 

In this case we have the following important values of R{r) ([3]), 

hm R(r) = 0, lim R(r) = oo. (4) 

Here, the solution represents space-times with a physical naked timelike singularity 
located at i? = j7] . It is easy to see that this singularity is physical because, the Ricci 
scalar R computed from the line element ([T]), 



-2 



diverges there. It is also easy to see that this singularity is naked with the aid of the 
quantity Q = 9°'^^^-§^- The roots of Q determine the presence of event horizons in 
spherical symmetric spacetimes [12]. For the Wyman solution Q is given by. 



g^(i-fl)-\i-^r. (6, 

The only root of Q is located at r = M + 7] which confirms that the singularity located 
at i? = is naked for all spacetimes in the present case. 

This singularity is sometimes called a 'central singularity' and is similar to that 
appearing in the 'extreme' Reissner-Nordstr0m black hole and in the negative mass 
Schwarzschild spacetime. From equation ([2]), the scalar field vanishes asymptotically as 
i? — > cxD and diverges at the singularity. 

Case M = 1] or X = 1. 

In this case we have R{r) = r and cr = 0. This last condition implies that the scalar 
field ([2]) vanishes and one gets the Schwarzchild solution. Here, the minimum of R{r), 
located at i? = 2M, is outside the domain of r and the line element ([T]) describes only 
the spacetime exterior to the event horizon. 

Case M > rj or 1 < X < oo. 

In this case we have the following important values of R{r) ([3]), 

limi?(r) = oo, /2^,„ = (M + r/)(-— -^)^' " \ lim i?(r) = oo.(7) 

Due to the fact that > — as stated just before ([2]), we have that in the present 
case 77 > 0. The case 77 = is well know in the literature as the Yilmaz- Rosen space-time 
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[TB] . In this case M > t], the physical radius R is never zero. If one starts with a large 
value of R, for large values of r, and starts diminishing R, reducing the values of r, 
one will reach the minimum value of R {Rmm ©) for r = M + t]. Then, R starts to 
increase again when we let r goes to 2r] until it diverges when r = 2t]. Therefore, we may 
interpret these spacetimes as wormholes connecting two asymptotically flat regions such 
that they have a minimum throat radius given by Rmm PJLJ- The spatial infinity 
(i? oo) of each asymptotically fiat region is obtained, respectively, by the limits: 
r — s> oo and r 2?]+. An important property of this space-time is that the scalar field 
([2]) is imaginary. The imaginary scalar field also known as ghost Klein- Gordon field [17] 
is an example of the type of matter called exotic by some authors [18]. It violates most 
of the energy conditions and is repulsive. This property helps explaining the reason why 
the collapsing scalar field never reaches R = 0. 

As mentioned above, in the rest of the paper we shall restrict our attention to the 
spacetimes representing naked singularities and wormholes. 



3. Timelike Geodesies 

3.1. Effective Potential 

We have four geodesic equations, one for each coordinate [1], 

d'^x" dx>^ dx^ 

where a = 0, 1, 2, 3, and x" represents, respectively, each of the coordinates {t, r, 9, 0). 
r is the proper time of the massive particle which trajectory is described by ([8]). 

The geodesic equation in the coordinate 9 tell us that, like in the Schwarzschild 
case, the geodesies are independent of 9, therefore we choose the equatorial plane to 
describe the particle motion {9 = tt/2). In the equatorial plane, the geodesic equation 
for may be once integrate to give. 



(l-^)''"0 = i?^0 = ^, (9) 



where the dot means derivative with respect to r and L is the integration constant that 
may be interpreted as the particle angular momentum per unit rest mass. This result 
means that is cyclic and its conjugated momentum (p^ = 0) is a conserved quantity. 
Also, may be written as a function of r. Likewise, in the equatorial plane, the geodesic 
equation for t may be once integrate to give, 

{l-^)H = E, (10) 

where E is the integration constant that may be interpreted as the particle energy per 
unit rest mass. This result means that t is cyclic and its conjugated momentum {pt = i) 
is a conserved quantity. Also, t may be written as a function of r. Instead of using the 
fourth geodesic equation for the r coordinate, we use the equation derived directly from 
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the line element ([T]), ds'^/dr'^ = —1 [IH]. There, we introduce the expressions of Q 
and t ffTOj) . in order to obtain the following equation which depends only on r, 

2 

+ V{ry = E\ (11) 




Where 



V{r) 



1-^1^ 



(12) 



^2 ( 

and V^(r) is the effective potential for the motion of massive particles in the Wyman 
geometry. The geodesic equation for the r coordinate plays the role of a control equation, 
where we substitute the solutions found with the other four equations in order to verify 
their correctness. 

In order to have a qualitative idea on the orbits of massive particles moving under 
the action of V^(r)^ it is important to draw it as a function of r. We may do it by, initially, 
computing the extremes of K(r)^. First, we calculate the first derivative of V{r)^ (1121) 
and find the roots of the resulting equation. That equation may be simplified by the 
introduction of the auxiliary quantities: x = (1 — 2?7/r), < x < 1; A = rf /L'^, 
0<A<oo;5=(A- l/2)/(A + 1/2), -I < B <l; C = 2A/(A + 1/2), < C < 2. 
Where A was defined before and the domains of A, B and C where determined by the 
fact that A is positive. The equation dV{r)/dr = 0, in terms of x, A, A, B and C is 
given by, 

ACx^-^ + B - Cx + x"^ = 0. (13) 

In order to identify the presence and nature of the roots of ( |T3l) . let us define the following 
two auxiliary functions, 

p(^x) = -x'^ + Cx - B , h{x)=ACx^-^. (14) 

Now, the values of x where the two curves p{x) and h{x) meet will be the roots of 
f fT3|) . p{x) is a set of parabolas which vertices are all located above the x-axis and 
with two roots. The larger one is located at x = 1 and the smaller one is located at 
X = (A — l/2)/(A + 1/2). For A > 1/2, the smaller root is positive; for A = 1/2, it is zero 
and for A < 1/2, it is negative. The precise nature of h{x) will depend on the value of 
A, present in the exponent of x, but whatever A one chooses, all values of h{x) will be 
located above the x-axis. For A > 1, h{x) diverges to -|-cxo when x — *• and goes to zero 
when X +00. For A < 1, h{x) goes to zero when x — > and diverges to +oo when 

X — s> +00. 

An important root of (fT3ll is defined by the value of x, say xq, where the two curves 
p{x) and h{x) f|T^ just touch each other, xq is an inflection point of V{rY f|T2l) . because 
there the second derivative of V{ry is also zero. Associated to Xq there is a value of the 
angular momentum L, say Lq, which originates an unstable particle orbit. Due to the 
fact that L is present, only, in the denominator of A in the expression of h{x) (IT^ . if 
one increases L, h{x) will assume smaller values for the same values of x. p{x) will not 
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be altered because L is not present in its expression. Therefore, Lq is the value of L for 
which h{x) just touches p{x), if one takes values of L greater than Lq, h{x) will start 
intercepting p{x) in two or more points. They will be extremes of maximum, 
minimum or inflection points. It is possible to compute the value of Lq in terms of Xq 
and the parameters A and rj. In order to do that, we consider, initially, the fact that the 
first derivatives of h{x) and p{x) (HM . in x = Xq, are equal and express ACxq'^ in terms 
of other quantities. Then, we use the fact that Xq is a root of f|T3|) and re-write that 
equation for x = xq and substitute there the value of ACxq^ just obtained. It gives rise 
to the following second degree polynomial equation in xq, 

where A 7^ 1. It has the following roots. 



2A^ + v/5A^ __ 2A^-v/5A^ 

2(A + 1/2)(A + 1)' 2(A + 1/2)(A + 1)' ^ ' 

where Xq > Xq . Due to the fact that we have two distinct values of Xq, Xq and Xq given 
by (ITBl) . we shall have, also, two distinct values of Lq, say Lo+ and Lq_. In order to 
obtain them, we introduce Xq and Xq , separately, in the equation that equates the first 
derivatives of h{x) and p{x) and express Lq in term of other quantities. It results in, 

r2 V'Kl - A) .2 ^^A(l-A) 

°+ (a;+)A[A-x+(A + l/2)]' ' {xq )'[X - Xq{X + 1/2)]' ^ ^ 

where Ll_ > Lq^. 

By definition A has the domain [0, cxd), but if one introduces values of A, over all 
this domain, in the expressions of Xq, Xq f fT6l) and Lq_,_, Lq_ f fTTl) . one will find some 
unphysical results. This means that one will have to impose some restrictions on the 
domain of A. These restrictions will have to take in account the following two conditions: 
(i) Xq and Xq varies in the range [0, 1], due to the definition of x; (ii) Lq_,_ and Lq_ must 
be positive. These conditions lead to the following distinct domains of A, depending on 
the extreme Xq or Xq one is using, 

x+ : ^ > (19) 

Therefore, this result tell us that V^(r)^ ( IT2|) behaves differently, depending on the 
value of A. We have the following three different regions: 

(i) A > 1/2. 

There will be just one inflection point located at Xq for Lq_,_. If one chooses values 
of > Lq^, one will find other extremes of V{r)'^ (1121) . which will be maximum or 
minimum points. 

(ii) 1/^5 < A < 1/2. 
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There may be two different inffection points. The first located at x'^ for = Lq_^. 
If one chooses values of > L^^, one will find other extremes of V"(r)^ f|T2l) . which 
will be maximum or minimum points. When one reaches = Lq_, one finds the other 
possible inflection point located at Xq . If one chooses values of > Lq_, one will 
find just one extreme of V{r)^ (fT2l) . which will be a minimum point. Even for the case 
< Lq_^^, there will be a minimum point. In fact, this minimum point will always be 
present for any value of L. Its presence can be understood because, here, p{x) has a 
negative root and h{x) is crescent and starts from x = 0. Therefore, these two curves 
will always intercept each other. 

(iii) A < 1/V5. 

There will be no inflection points but there will be always a minimum point. The 
presence of this minimum point can be understood in the same way as the one in the 
previous case. 

As we saw in the previous section. Sec. |2l naked singularities and wormholes are 
characterized by certain subdomains of A. Therefore, based on the above result, we 
may draw the effective potentials for naked singularities and wormholes. From these 
effective potentials, we shall be able to describe qualitatively the different orbits of 
massive particles. 

3.2. Effective Potential for Naked Singularities 

As we saw, in Sec. [21 the naked singularities are obtained for < A < 1. We may still 
divide the naked singularities in two classes, due to the behavior of V(r)^ (fT2l) when 
r —>■ 2?]+. For < A < 1/2, V"(r)^ ( fT2l) assumes the following values when r —>■ 2?7+ and 
r +00, 

lim V(rf = oo, lim V(rf = 1. (20) 

Due to the fact that these limits are consistent with an asymptotically flat naked 
singularity located at r = 27], we call this class of ordinary naked singularities. On 
the other hand, for 1/2 < A < 1 the limit of V"(r)^ ( |T2|) when r — > 2?7+ is zero. Due to 
fact that this result is not consistent with a naked singularity located at r — 27], we call 
this class of anomalous naked singularities. Then, in what follows we shall restrict our 
attention to the class of ordinary naked singularities with < A < 1/2. It is important 
to mention that observing the scalar field expression ([2]), which in this case is real, one 
can see that it diverges to — oo as r — 2?7+. 

Taking in account the results of Subsec. 13.11 (ii) and (iii) the effective potential 
V{r)'^ (|T2l) may have several different shapes depending on the value of L^. Here, the 
inflection points will be located at Xq and Xq f|T6|) and the relevant angular momenta 
will be Ll_ and L^^ 1^7^). 

• For < -^0+! y{f^y ( fT2l) has one minimum point. In terms of x it is located in the 
range {0,Xq ) or in terms of r it is located in the range {2r] , Tq = 2r]/ (1 — Xg )). For 
E"^ < 1, the massive particles orbit around the naked singularity. If the massive 
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particle is located exactly at the minimum point the orbit is circular and stable. 
For ^2 > 1, the massive particles come in from infinity reach the infinity potential 
wall near the naked singularity and return to infinity without ever reach the naked 
singularity. 

• For = Lq^, V{ry f|T2|) has two extreme points, a minimum located in the range 
(2?7 , r^") and an infiection point located at r = Tq = 21] / {1 — Xq) . In this point the 
massive particles have unstable circular orbits. The other possible trajectories for 
the massive particles are exactly as in the previous case. 

• For Lq_,_ < < Lq_, V{ry f|T2l) has three extreme points, a minimum located 
in the range {2ri , r^), a maximum located in the range {r^ , r^) and another 
maximum located in the range {tq , oo). If the massive particles are exactly in the 
maximum points they have unstable circular orbits. The other possible trajectories 
for the massive particles are exactly as in the first case. 

• For = -Lq_, K(r)^ ( |T2l) has two extreme points, a minimum located in the 
range {tq , oo) and an infiection point located at r = Tq . In this point the 
massive particles have unstable circular orbits. The other possible trajectories for 
the massive particles are exactly as in the first case. 

• For > Lq_, V{ry (1121) has one minimum point located in the range {tq , oo). 
The possible trajectories for the massive particles are exactly as in the first case. 

It is important to mention that, for sufficiently large values of E, any massive 
particles will come in from infinity and it will get refiected by the infinity potential wall 
near the naked singularity. Then, it will return to infinity without ever reach the naked 
singularity. This result is similar to the one found in [3] for timelike geodesies of the 
naked singularity present in the Reissner-Nordstrom spacetime. Let us consider, as an 
example, the case where M = 10 and A = 0.45. Therefore, we have Lq_,_ = 954.9308516 
and Xq = 0.1875874406 which in terms of r is given by Vq = 54.70674220. We also 
have, Lq_ = 960.0309665 and Xq = 0.1064234487 which in terms of r is given by 
rQ = 49.73770224. In the Figure 1, we plot five different effective potential diagrams 
from V{ry ( fT2l) . one for each of the five cases discussed above. They have the 
five different values of L^: 951, 954.9308516, 957, 960.0309665 and 962. The naked 
singularity is located at r = 44.44444444. 



3.3. Effective Potential for Wormholes 

As we saw, in Sec. [2l the wormholes are obtained for 1 < A < oo. In the spatial 
infinity of each asymptotically fiat region the effective potential V"(r)^ ( 1121) assumes the 
following values, 

lim V(rf = 0, lim V(rf = 1. (21) 

Although the above limits give consistent values for the effective potential at the two 
spatial infinities, the first limit is not consistent with the value of the Ricci scalar 
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Figure 1. Five different effective potential diagrams from V{r)'^ (US), for massive 
particles moving in a naked singularity with M — 10 and A = 0.45. They have the 
five different values of L'^: 951, 954.9308516, 957, 960.0309665 and 962. The naked 
singularity is located at r = 44.44444444. 



evaluated at r = 27], for wormholes with 1 < A < 2. For this class of wormholes the 
Ricci scalar R ([5]), diverges to oo when we take the limit r 21]^ or gives a positive 
constant when A = 2. These values are not consistent with the first limit of l^(r)^ in 
fl2T]) and the idea of an asymptotically flat spatial region. Therefore we call this class of 
wormholes of anomalous wormholes. For the rest of the wormholes where 2 < A < oo, 
the limit of R ([5]) when r — >■ 2r]^ is zero. This value is consistent with the first limit of 
V{ry in (EH) and the idea of an asymptotically fiat spatial region. Then, in what follows 
we shall restrict our attention to that class of ordinary wormholes with 2 < A < oo. It 
is important to mention that observing the scalar field expression ([2]), which in this case 
is imaginary, one can see that it diverges to — oo as r ^ 2?]+. Since there is no physical 
singularity at r ^ 2ri^ for the ordinary class or wormholes, this result indicates that 
the coordinates we are using might not be suitable to describe this asymptotically flat 
region. 

Taking in account the results of Subsec. 13.11 (i) the effective potential V"(r)^ ( fT2l) 
may have three different shapes depending on the value of L^. Here, the inflection point 
will be located at Xq (fT6|) and the relevant angular momentum will be L^^ (ITTl) . 

• For < -^0+) V{ry (fT2l) has no extreme points. In this case, there is no stable 
orbits. For sufficiently high energies the massive particles may travel from one 
asymptotically flat region to the other. In fact, this type of orbit is also present in 
the other two cases considered below. 

• For = Lq_^, V{ry f|T2l) has one inflection point, located at Xq (fT6l) or in terms 
of r, = 2rj/{l — Xq). In this point the massive particles have unstable circular 
orbits. 
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• For > Lq_^, V{ry (fT^ has two extreme points, a maximum located at 
2?7 < r < Tq" and a minimum located at < r. There are closed and open 
orbits depending on the values of the total energy and angular momentum of the 
massive particles. 

Let us consider, as an example, the case where M = 1 and A = vTOOO. Therefore, 
we have Lq_^ = 12.41266 and Xq = 0.98799 which in terms of r is given by Tq = 5.26747. 
In the Figure 2, we plot three different effective potential diagrams from y(r)^ ( IT2l) . 
one for each of the three cases discussed above. They have the three different values of 
L^: 10, 12.41266 and 14.4. The spatial infinities of each asymptotically flat region are 
located at r = 0.06325 and r — > oo. 
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Figure 2. Three different effective potential diagrams from V{r)'^ (fT2| . for massive 

particles moving in a wormhole with M — 1 and A = VIOOO. They have the 

three different values of I?: 10, 12.41266 and 14.4. The spatial infinities of each 
asymptotically flat region are located at r = 0.06325 and r ^ co. 

It is important to mention, at this point, that the same study of the qualitative 
features of the orbits of massive particles in Wyman geometry in terms of the proper 
time r, as we just did, can also be done in terms of the coordinate time t. In order to 
do it, we re- write equation (fTTj) in the following way, 

[e^-^ ^V{rf = E\ (22) 

where we used the relationship, dt/dr = E/{1 — 2ri/r)^, and we introduced the 
generalization of the Wheeler's tortoise coordinate r* pQ, given by. 
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Since, V{r) and E in (1221) are the same as the ones in fill I) , the turning points of V{r), 
for a given value or E, are the same whether we use t or r. On the other hand, the 
geodesies may have different properties depending on the time variable we are using. 
We shall see some similarities and differences in the next section. 

3.4- Radial timelike geodesies 

3.4-1 ■ Description in terms of r Unfortunately, there is not an analytic expression for 
the general timelike geodesies given by the solutions of Eqs. fl MTTl) . even the numerical 
solutions are very complicated. On the other hand, one may restrict his attention to 
the case of radial timelike geodesies, where the massive particle moves only along the 
radial and time directions. It means that 6 = (p = L = and (ITTll reduces to. 



Although (12^ is much simpler than ( fTTl) . one still cannot integrate it to find an 
analytic expression of r as a function of r. The best one can do is finding a numerical 
solution for the integral and plot that solution in a graph of r versus r. We did that 
and we shall present the results below, for the naked singularities and wormholes. 

Naked singularities In this case, we integrated (12^ for many different values of E, M 
and rj. We chose values of M and rj compatibles with naked singularities. We found 
that, for E > 1, the geodesies are all well behaved when r 2ri. In fact, r goes to zero 
in that limit. The geodesies are such that when r is large, r tends to a linear function 
of the proper time r. For E < 1, the geodesies are also well behaved when r ^ 2r) 
and r also goes to zero in that limit. On the other hand, they do not extend to large 
values of r. It is clear from that the massive particle is subjected to a potential 
of the form: (1 — 2r]/r)^'^^^. This potential, increases from zero at r = 2r] and tends to 
one when r oo. Therefore, massive particles with a total energy E < 1, get reflected 
by the potential. The value of r, where the particle gets reflected by the potential is 
obtained by solving the equation: E'^ = (l — 2r]/r)^/^. Then, based on the above results, 
we may conclude that it always takes a finite proper time interval to travel from any 
finite value of r to the singularity located at r = 2r]. We notice, also, that for radial 
timelike geodesies there is not an infinity potential wall near the naked singularity and 
the massive particles can reach it. In the Figure 3, we show as an example a geodesic for 
the case: E = 1, M = 1 and t] = 3. The naked singularity is located at r = 6 and the 
geodesies extends to large values of r. In the Figure 4, we show as an example a geodesic 
for the case: E = 0.5, M = 1 and rj = 3. The naked singularity is located at r = 6 and 
the geodesies gets reflected at r = 6.095238095. One cannot see the geodesic returning 
to r = 6, from this figure, because we have considered only the outgoing geodesic. 

Wormholes In this case, we integrated for many different values of E, M and ry. 
We chose values of M and rj compatibles with wormholes. We found that, for E > 1, the 




(24) 
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Figure 3. Radial geodesic, described in terms of r, for a massive particle moving in 
a naked singularity with E = 1, M = 1 and = 3. The naked singularity is located at 
r = 6. 
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Figure 4. Radial geodesic, described in terms of t, for a massive particle moving in 
a naked singularity with E = 0.5, M = 1 and r/ = 3. The naked singularity is located 
at r = 6. 

geodesies are all well behaved when r — > 21]. In fact, r goes to zero in that limit. The 
geodesies are such that when r is large, r tends to a linear function of the proper time 
r. For E < 1, the geodesies are also well behaved when r — > 2?] and r also goes to zero 
in that limit. On the other hand, they do not extend to large values of r because they 
get reflected by the potential, as in the naked singularity case. The value of r, where 
the particle gets reflected by the potential is obtained in the same way we proceeded in 
the naked singularity case. Here, as in the naked singularity case, for all cases studied 
it always takes a finite proper time interval to travel from any finite value of r to the 
spatial infinity located at r = 21]. In the Figure 5, we show as an example a geodesic 
for the case: E = 1, M = 3 and i] = 1. The spatial infinities of each asymptotically fiat 
regions are located at r = 2 and r — > cxo and the geodesies extends to large values of r. 

3.4-2. Description in terms oft Here, as in the case where the geodesies are described 
in terms of r, the treatment of the general situation is very complicated. Therefore, we 
shall, also, restrict our attention to the particular situation of radial timelike geodesies. 
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Figure 5. Radial geodesic, described in terms of t, for a massive particle moving in a 
wormhole with E ^ I, M = 3 and 77 = 1. The spatial infinities of each asymptotically 
flat regions are located at r = 2 and r 00. 



From Eqs. f|T2l) . fl22|) and fl23|) . appropriately written for the radial case, we may obtain 
the expression which, after integration, gives r as a function of t. 



Even with this simplification, as in the case with r, we cannot integrate ( !25l) to find an 
analytic expression of t as a function of r. The best one can do is finding a numerical 
solution for the integral. We did that for both cases of the naked singularities and 
wormholes. 

Naked singularities In this case, we integrated (l25l) for many different values of E, M 
and rj. We chose values of M and rj compatibles with naked singularities. Qualitatively, 
the geodesies behave very much like the description in terms of r. For E > 1 and E < 1, 
t goes to zero when r ^ 2r]. For E > 1, when r is large, r tends to a linear function of 
t. For E < 1, the geodesies do not extend to large values of r. As we have mentioned 
above, since the potential is the same as in the description in terms of r, the turning 
points are also the same. Nevertheless, there is a quantitative difference between the two 
descriptions. We observe that, although, in both descriptions r is an increasing function 
of the time variable, it grows slower with t than with r. The specific difference in the 
rate of growth between the two descriptions depends on the values of the parameters M 
and rj. Here, we may, also, conclude that it always takes a finite time coordinate interval 
to travel from any finite value of r to the singularity located at r = 2ri and there is not 
an infinity potential wall near the naked singularity and the massive particles can reach 
it. 

Wormholes In this case, we integrated (125|) for many different values of E, M and 77. 
We chose values of M and t] compatibles with wormholes. In analogy with the case 
where the geodesies are described by r, for E > 1, when r is large, r tends to a linear 




(25) 
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function of t. On the other hand, for E < 1 they do not extend to large values of r. 
As we have mentioned above, since the potential is the same as in the description in 
terms of r, the turning points are also the same. We found, here, some differences with 
respect to the description in terms of r. The most important is that, for any value of 
E, the geodesies are not well behaved when r — > 2?]+. In fact, t goes to — oo in that 
hmit. As another difference, we observe that, although, in both descriptions r is an 
increasing function of the time variable, it grows slower with t than with r. The specific 
difference in the rate of growth between the two descriptions depends on the values 
of the parameters M and r). Due to the divergence in t when r — > 2r)+, we conclude, 
here, that it always takes an infinite time coordinate interval to travel from any finite 
value of r to the spatial infinity located at r = 21]. In the Figure 6, we show as an 
example a geodesic for the case: E = 1, M = 3 and rj = 1. The spatial infinities of each 
asymptotically fiat regions are located at r = 2 and r —>■ oo and the geodesies extends 
to large values of r. 



4. Null Geodesies 

4.1- Effective Potential 

The null geodesies for the Wyman solution are derived almost in the same way the 
timelike geodesies were derived in the previous section. The only difference is that, 
here, the null line clement contributes a different additional equation to the four geodesic 
equations. The new equation reads: ds'^/dx^ = 0, where x is the affinc parameter used 
in the present case. Therefore, proceeding exactly as in the previous section we obtain 
the following effective potential equation. 




-so - 



Figure 6. Radial geodesic, described in terms of t, for a massive particle moving in a 
wormhole with E = 1, M = 3 and r] = 1. The spatial infinities of each asymptotically 
flat regions are located at r = 2 and r — > oo. 




(26) 
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where 

V{rf = r-\l-'^)'^~' (27) 

and b = L/E is the photon impact parameter. V"(r)^ fl271) is the effective potential for 
the motion of the photon in Wyman geometry. 

V^(r)^ (I27|) has only one extreme point at r = 2M + r]. Due to the domain of r it can 
only exists if 2M > rj or X > 1/2 and when it exists it is a maximum point. Therefore, 
we have three different cases: (i) A > 1/2, the effective potential has a maximum point 
at r = 2M + rj and goes to zero at r = 2ri;{n) A = 1/2, the effective potential has no 
extreme points and diverges to oo as r ^ 0; (iii) A < 1/2, the effective potential has no 
extreme points and diverges to oo as r — 2?]. 

As we saw in Sec. |2l naked singularities and wormholes are characterized by certain 
subdomains of A. Therefore, based on the above result, we may draw the effective 
potentials for naked singularities and wormholes. From these effective potentials, we 
shall be able to describe qualitatively the different orbits of photons. 



4-2. Effective Potential for Naked Singularities 

As it was mentioned in Subsec. 13.21 we are only concerned, here, with the class of 
ordinary naked singularities characterized < A < 1/2. For this class, V"(r)^ ( |T2l) 
assumes the following values when r — > 2?7+ and r +oo, 

lim V(r? = oo, lim V(r? = 0. (28) 

\^(r)^ fl271) has no extreme points. In this case, whatever the impact parameter b the 
photons come in from infinity get reflected by the infinity potential wall near the naked 
singularity and return to infinity without ever reach the naked singularity. 

Let us consider, as an example, the case where M = 0.3, rj = 1 and A = 0.3. In the 
Figure 7, we plot the effective potential diagram from V^(r)^ 0271) . for these parameter 
values. The naked singularity is located at r = 2. 



4-3. Effective Potential for Wormholes 

As it was mentioned in Subsec. 13.31 we are only concerned, here, with the class of 
ordinary wormholes characterized by 2 < A < cxo. In the spatial infinity of each 
asymptotically flat region the effective potential V{r)^ (1271) assumes the following values, 

lim V(r? = 0, lim V(r? = 0. (29) 

V{ry (p7|) has a maximum point located at r = 2M + rj. In this point the photons 
have unstable circular orbits. For sufficiently small impact parameter b the photons may 
travel from one asymptotically flat region to the other. Otherwise, they come in from 
spatial infinity of an asymptotically flat region get reflected by the effective potential 
and return to spatial infinity in the same asymptotically flat region. 
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Figure 7. Effective potential diagram from V{r)'^ (l?7)) . for photons moving in a naked 
singularity with M = 0.3, 77 = 1 and A = 0.3. The naked singularity is located at r = 2. 

Let US consider, as an example, the case where M = 30, 77 = 1 and A = 30. In the 
Figure 8, we plot the effective potential diagram from ^(r)^ flTTI) . for these parameter 
values. The spatial infinities of each asymptotically flat regions are located at r = 2 
and r — > 00. 



0,00003- 



0,00002- 

0,00001 - 



0,00000 

50 100 150 200 250 300 



Figure 8. Effective potential diagram from V{r)'^ (f27| . for photons moving in a 
wormhole with M = 30, ry = 1 and A = 30. The spatial infinities of each asymptotically 
flat regions are located at r = 2 and r 00. 



4-. 4- Radial null geodesies 

4.4-1 ■ Description in terms of x The description of the null geodesies in terms of the 
affine parameter x is trivial, when we consider radial geodesies {9 = (j) = L = {]). That 
is the case because the effective potential equation for the radial motion of photons 
reduces to (dr/dx)^ = E"^- It means that the effective potential is zero and the photons 
move freely following straight lines: r = ±Ex + q, where q is an integration constant 
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and the + and — signs mean, respectively, outgoing and ingoing geodesies. This result is 
valid for photons moving in the naked singularity as well as in the wormhole geometries. 



4-4-2- Description in terms oft The description of the null geodesies in terms of the 
time coordinate t may be derived with the aid of ( ITOl) . appropriately written in terms 
of X, and the equation {dr/dxY = Then, we have: (dt/dr)'^ = [dt / dxY I {dr / dxY = 
(1 — 2r] / r)^~'^^\ From this equation, we obtain t as a function of r after performing the 
following integral. 



t = ± 



2r] 

r 



dr 



Fortunately, we can analytically integrate (!30|) to find. 



t = ± 



, /Tn 2M(M + r]) /M 2r?N 
r + 2Mln(-) ^ ^ g^J— + 2 , 1 , 1 ; 2 , 3 ; ^) 



(30) 



-Q.(31) 



where the + and — signs mean, respectively, outgoing and ingoing geodesies. Q is an 
integration constant and F is an hypergeometric function. If we introduce the condition 
M = in ( l3Ti) . we obtain: t = r + 2Mln(r — 2M) + Ql, where Q\ is an integration 
constant. That expression is the well-known radial, outgoing null geodesies for the 
Schwarzschild geometry parametrized by the time coordinate [T]. 

For large values of r (|3T1) indicates that r becomes a linear function of t, for both 
naked singularities and wormholes. On the other hand, a great difference appears 
between naked singularities and wormholes in the limit when r — > 27]+. Observing 
(!3T|) . we see that the limit of t when r 2ri^ depends crucially on the behaviour of 
3F2 when r — > 2//+. From the theorem 2.1.2, page 62 in [20], we learn that in the limit 
r 2174., 3F2(M/r7 + 2, 1, 1; 2, 3; 2?7/r) converges absolutely if Mjr] < 1 and diverges 
otherwise [M/t] > 1). Therefore, for naked singularities t will converge for a finite value 
and for wormholes it will diverge to —00. It is actually possible to compute the limit of 
3^2(71^/77 + 2, 1, 1; 2, 3; 2ri/r) when r — > 2ri^. In order to do that, we must use: equation 
(2.2.2), page 67 in [20j, 2Fi{l ,1 -,2 ; z) = -i In | 1 - z | and r(a;)r(l - x) = 7r/sin(7rx). 
After an integration we find the following result. 



3F2(A + 2, 1, 1; 2, 3;1) 



2ri 




(32) 



M + r] 

where ip{x) = r'{x)/r{x) is the digamma function, T{x) is the gamma function and 
7 = 0.5772156649 ... is the Euler constant. Therefore, introducing this result in (13T!) 
we obtain the value of t when r 2ri+, 



lim t 

r— »2»7+ 



2{ri -M)+2M 



+ 7 



+ Q 



(33) 



Based on the above results, we may conclude that for the case of naked singularities 
it always takes a finite time coordinate interval to travel from any finite value of r to 
the singularity located at r = 2//. On the other hand, it always takes an infinite time 
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coordinate interval to travel from any finite value of r to the spatial infinity located 
at r = 2?7, in the case of wormholes. Both results are qualitatively analogous to the 
corresponding ones found for the radial timelike geodesies. Finally, we notice, also, that 
for null radial geodesies there is not an infinity potential wall near the naked singularity 
and the photons can reach it. 

5. Conclusions. 

In the present work, we studied qualitative features of orbits of massive particles and 
photons moving in the naked singularity and wormhole spacetimes of the Wyman 
solution. These orbits are the timelike geodesies for massive particles and null geodesies 
for photons. Combining the four geodesic equations with an additional equation derived 
from the line element, we obtained an effective potential for the massive particles and 
a different effective potential for the photons. We investigated all possible types of 
orbits, for massive particles and photons, by studying the appropriate effective potential. 
We noticed that for certain values of cr^ > 0, there is an infinity potential wall that 
prevents both massive particles and photons ever to reach the naked singularity. This 
result is similar to the one found in [3] for timelike geodesies of the naked singularity 
present in the Reissner-Nordstrom spacetime. We noticed, also, that for certain values 
of — Mq < cr^ < 0, the potential is finite everywhere, which allows massive particles 
and photons moving from one wormhole asymptotically flat region to the other. We 
also computed the radial timelike and null geodesies for massive particles and photons, 
respectively, moving in the naked singularities and wormholes spacetimes. 
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